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PREFACE 


This  filial  report  is  a description  of  the  technology  used  and 
results  obtained  in  a test  intended  to  verify  the  feasibility  of  using 
laser  Doppler  velocimetry  for  measurement  of  low  speed  winds  in 
fog.  The  test  was  conducted  by  personnel  of  Lockheed  Missiles  & 
Space  Company,  Inc.,  Huntsville  Research  & Engineering  Center. 
Lockheed-Huntsville  personnel  contributing  to  this  effort  were  E.  J. 
Gorzynski,  G.M.  Miller  and  K.R.  Shrider.  The  dedicated  efforts 
during  the  test  phase  of  Dr.  Alan  I.  Weinstein  and  Frederick  J. 
Brousaides  of  the  Air  Force  Geophysics  Laboratory  are  gratefully 
acknowledged. 
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Section  1 

INTRODUCTION  AND  SUMMARY 


1.1  BACKGROUND 

Significant  effort  is  currently  being  devoted  to  development  of  instru- 
mentation to  remotely  sense  atmospheric  flow  phenomena.  Some  of  the 
avenues  being  pursued  are  active  and  passive  acoustic  sensors,  optical 
sensors,  and  radio  methods.  A useful  survey  of  such  methods  is  presented 
in  Ref.  1.  Two  advantages  of  remote  sensors  are  that  flow  conditions  can  be 
ascertained  in  regions  of  space  where  it  would  not  be  convenient  to  locate 
instrumentation  hardware,  and  no  interference  with  the  flow  at  the  point  of 
interest  is  introduced  by  their  use.  The  laser  Doppler  velocimeter  (LDV) 
is  a particularly  attractive  device  for  remote  sensing  of  atmospheric  phe- 
nomena. In  the  LDV  system,  the  laser  radiation  backscattered  by  moving 
particulates  in  the  atmosphere  is  used  to  determine  the  velocity  of  the  flow. 
Since  it  is  possible  to  direct  the  laser  focal  volume  at  a selected  sequence 
of  points  in  space,  data  from  a scanning  LDV  system  cam  be  used  to  deter- 
mine the  velocity  field  rapidly  and  over  a range  of  altitudes.  A CC>2  laser 
Doppler  velocimeter  system  has  the  following  advantages  over  other  remote 
sensing  techniques:  (1)  the  position  of  the  volume  within  which  velocity  is 
sensed  can  be  varied  with  ease  as  only  optic  pointing  and  focusing  operations 
are  involved;  (2)  the  ambient  aerosol  provides  a sufficient  scattering  target; 
and  (3)  the  sensing  mechanism  is  non -mechanical  which  results  in  the  potential 
for  a high  frequency  turbulence  sensor. 

The  feasibility  of  using  a laser  Doppler  velocimeter  (LDV)  system  for 
the  remote  sensing  of  low  altitude  winds  and  for  the  detection  and  tracking  of 
aircraft  wake  vortices  has  been  demonstrated  (Refs.  2 through  6).  The  par- 
ticular application  for  which  the  test  described  herein  was  conducted  is  the 
measurement  of  wind  in  fog.  The  accurate  measurement  of  wind  in  fog  is 
necessary  for  successful  application  of  fog  dispersal  systems  to  dissipate 
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fog  at  military  bases.  The  accurate  measurement  of  low  speed  wind  is 
particularly  important. 

t 

1.2  PROGRAM  OBJECTIVES 

The  primary  objective  of  the  test  was  to  verify  the  ability  of  the  LDV  « 

to  measure  low  speed  winds  in  both  advection  fog  and  radiation  fog.  It  was 
also  desired  to  define  particular  potential  system  or  operating  modifications 
which  could  be  employed  to  improve  low  speed  wind  measurement  capability 
in  a system  specifically  designed  for  such  purposes.  An  adjunct  objective 
was  an  evaluation  of  the  feasibility  of  measuring  visibility  with  the  LDV  by 
measuring  the  attenuation  of  the  laser  signal.  The  results  of  that  feasibility 
evaluation  are  not  reported  herein.  However,  the  adjunct  purpose  of  the  test 
is  discussed  because  it  affected  the  manner  in  which  the  wind  measurement 
data  were  taken. 

1.3  REPORT  FORMAT 

A brief  description  of  laser  Doppler  velocimetry  and  the  Lockheed- 
Huntsville  Mobile  Atmospheric  Unit  is  presented  in  Section  2.  Section  3 
contains  a description  of  the  mathematical  algorithms  used  to  measure 
winds.  The  tests  conducted  at  the  Air  Force  Geophysics  Laboratory  Test 
Site  at  Otis  AFB,  Massachusetts  are  described  in  Section  4 with  a presenta- 
tion and  discussion  of  results. 
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Section  2 

LASER  DOPPLER  SYSTEM  DEVELOPMENT 


2.1  SYSTEM  DESCRIPTION 

An  LDV  wind/vortex  sensor  senses  air  movement  by  measurement  of 
the  Doppler  frequency  shift  of  laser  radiation  backscatter ed  by  the  atmos- 
pheric aerosol.  An  instrument  must  incorporate  means  to  transmit  the  laser 
adiation  to  the  region  of  interest,  collect  the  radiation  scattered  from  the 
atmospheric  aerosol  and  to  photomix  the  scattered  radiation  and  a portion  of 
the  transmitted  beam  on  a photodetector.  The  difference  between  the  trans- 
mitted frequency  and  the  returned  frequency  is  the  Doppler  shift  frequency. 
The  Doppler  frequency  shift  signal  is  generated  at  the  photodetector  and  is 
translatable  into  an  along -optic  axis  wind  velocity  component  using  appro- 
priate electronics.  The  magnitude  of  the  Doppler  shift,  Af,  is  given  by  the 
equation 


Af  = y | V[  cosQ 
A 


where 


( V|  = the  velocity  vector  in  the  region  being  sensed 

X = the  laser  radiation  wavelength,  and 

0 = the  angle  subtended  by  the  velocity  vector  and 

the  optic  system  line  of  sight. 

A Doppler  shift  of  188  MHz  results  per  m/sec  of  line -of- sig ht  velocity  com- 
ponent. Thus  measurement  of  the  Doppler  shift  frequency,  Af,  yields  directly 
the  line-of-sight  velocity  component  I vl  cos9.  Some  typical  advantages  of  the 
laser  Doppler  method  are:  (1)  the  Doppler  shift  is  a direct  absolute  measure 
of  the  velocity  (for  example,  the  hot  wire  yields  velocity  via  a cooling  effect 
on  the  wire),  (2)  the  ease  with  which  the  position  of  the  sensing  volume  can 
be  varied  (optics  pointing  and  focusing  operations  only  being  involved);  (3) 
the  ambient  aerosol  provides  sufficient  scattering,  thus  enabling  operation 
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in  "clear  air"  conditions;  and  (4)  the  ambient  aerosol  tracer  has  a small 
inertia  and  responds  quickly  to  variations  in  airspeed  and  is  thus  a good 
turbulence  indicator. 

A useful  instrument  must  also  incorporate  means  of  scanning  the  sys- 
tem's sensing  volume  in  the  desired  manner  and  to  also  effect  the  required 
signal  processing,  on-line  read-out  and  permanent  recording  requirements. 

The  hardware  implementation  of  the  field  laser  Doppler  unit  utilized 
during  this  investigation  is  discussed  in  the  following  subsections.  The 
overall  configuration  is  summarized  in  Fig.  2-1. 

2.1.1  Basic  LDV  Optical  System 

The  optical  system  is  of  a monostatic  design  and  utilizes  a continuous 
wave  laser.  This  arrangement  depends  on  focusing  the  transmitter  telescope 
at  the  location  of  interest  for  its  spatial  resolution  property.  Details  of  the 
optical  arrangement  are  shown  in  Fig.  2-2. 

Specifically,  a horizontally  polarized,  20 -watt,  continuous  wave  CO^ 
laser  beam  (10.6  micron  wavelength)  emerges  from  the  laser  (T)  and  is 
deflected  90  degrees  by  a mirror  (T).  The  approximately  6 mm  diameter 
beam  then  passes  through  a Brewster  window  (4)  and  a CdS  quarter  wave- 
plate  (?)  which  converts  it  to  circular  polarization.  The  beam  impinges  on 
the  secondary  mirror  (V)  and  is  expanded  and  reflected  into  the  primary  mir- 
ror (30  cm  diameter)  (V)and  then  focused  out  into  the  atmosphere.  A small 
portion  of  the  original  laser  beam,  is  reflected  by  the  secondary  mirror  and 
the  Brewster  window  (?)  and  is  used  as  a reference  frequency  on  the  photo- 
detector (lO)  . Energy  scattered  by  aerosols,  at  the  focal  volume  (?)  is 
collected  by  the  primary  mirror  (?),  collimated  by  the  secondary  (?) , and 
passed  through  the  quarter  waveplate  (?).  The  quarter  waveplate  changes 
the  polarization  of  the  aerosol  backscattered  radiation  from  circular  to 
vertical  linear  polarization.  The  vertically  polarized  beam  is  approximately 
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IMir  ror 
(Mirror 

(Brewster  Window 
'Quarter  Wave  Plate 
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Primary  Mirror 
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Photodetector 
Preamplifier 
Spectrum  Analyzer 


Fig.  2-2  - Optical  Component  Configuration  of  the  Lockheed  LDV 


78%  reflected  off  the  Brewster  window  (4).  After  passing  through  the  col- 
lecting lens  (?)  the  two  beams  are  photomixed  on  the  detector  (lO)  in  a 
heterodyne  configuration.  The  electrical  output  of  the  detector  (To)  is 
amplified  (0)  with  a 5 MHz  bandwidth,  20  dB  gain  low'  noise  type  preamp- 
lifier and  fed  into  a spectrum  analyzer  (T2) . 


2.1.2  Optic  Scanning  System 

In  order  to  provide  the  flexibility  required  to  operate  the  various  re- 
quired modes,  a scanning  arrangement  as  shown  in  Fig.  2-3  is  utilized.  The 
required  modes  of  operation  include  vortex  tracking  (not  required  for  the 
measurements  described  in  this  document)  and  velocity  azimuth  display  (VAD) 
for  measurement  of  atmospheric  wind.  The  mirror  assembly,  AB , can  be 
rotated  about  the  vertical  axis  for  scanning  in  azimuth  necessary  for  the  VAD 
(also  called  conical  scan  mode  of  operation).  Mirror  A is  adjusted  to  control 
the  elevation  angle  of  the  beam,  thus  controlling  the  cone  angle  of  the  conical 
scan.  The  scanning  hardware  as  deployed  on  the  mobile  van  is  shown  in 
Fig. 2-4. 

Range  scanning  of  the  system's  focal  volume  is  accomplished  by  vary- 
ing the  distance  between  the  telescope  secondary  mirror,  E,  and  the  primary 
mirror , D . Thi s is  effected  by  varying  the  position  of  the  mirror,  E,  in  a 
controlled  manner  by  an  electric  motor/optical  encoder  combination. 

The  operator  inputs  for  the  scanning  system  are  made  through  a control 
panel  incorporating  thumbswitch  controls  and  EED  monitors.  The  system's 
scan  capabilities  is  summarized  in  Table  2-1  and  Fig . 2-5, respectively . 

2.1.3  Signal  Processing  System 

The  Doppler  frequency  shift  of  the  photodetector  output  is  ascertained 
through  use  of  a sampled  spectrum  analyzer  which  provides  frequency  spectra 
(intensity  of  returned  signal  as  a function  of  Doppler  shiJt.)  at  a rate  of  <0  sig- 
natures per  second.  A typical  Doppler  wind  signature  is  shown  in  Fig.  2-6. 
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Table  2-1 


RANGE 


ELEVATION 


VAD  MODE 


MULTIMODE 


AC  C U RAC V 


SCAN  CAPABILITIES 


Maximum  Limit:  100  to  800  m 
Minimum  Limit:  16  to  650  m 
Scan  Frequency:  0.1  to  6.9  Hz 

Maximum  Limit:  10  to  90  deg 
Minimum  Limit:  0 to  90  deg 
Hard  Limit:  3 deg  with  Override 
Scan  Frequency:  0.1  to  0.5  Hz 

Measurement  Altitude:  10  to  640  m 
Measurement  Time/Altitude:  5 sec 
Sample  Rate:  1 to  7 Cycles 
Number  of  Altitudes:  8 

Elevation  Coverage:  3 to  90  deg  Upwind  and 
Downwind 

Scan  Plane  Azimuth:  360  deg 
Vertical  Line  Scan:  16  to  640  m 
Overhead  Arc  Scan:  90  cleg  Coverage 

Range:  0.5  m at  500  ft 
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To  yield  a line-of- sight  velocity  estimate,  a voltage  is  made  available  which 
has  the  same  time  behavior  as  the  Doppler  shift,  f^,  as  given  by  the  peak  of 
the  spectrum. 

The  implementation  of  this  technique  is,  in  essence,  a recursive  com- 
parison method.  The  spectrum  analyzer  scan  is  driven  by  a sawtooth  voltage 
derived  from  a D/A  converter,  the  input  to  which  is  counter  clocked  at  a con- 
stant rate,  hence  the  digital  number  output  of  the  counter  represents  frequency 
on  a linear  scale.  At  each  new  count,  the  spectrum  analyzer  output  is  con- 
verted to  a digital  representation  by  an  A/D  converter,  and  the  binary  number 
representing  the  current  sample  is  compared  with  that  obtained  on  the  pre- 
vious count.  If  the  current  one  is  the  larger  of  the  two,  it  is  saved  by  storing 
in  a latch,  along  with  the  binary  number  representing  its  frequency;  if  it  is 
smaller,  the  previous  one  is  retained  until  the  next  comparison.  This  process 
is  continued  for  the  entire  sweep.  It  is  evident  that  the  number  remaining  in 
the  frequency  store  latch,  when  the  sweep  is  completed,  corresponds  to  the 
highest  signal  power  observed,  i.e.,  the  peak  of  the  spectrum.  At  the  end  of 
each  sweep,  the  new  peak  frequency  found  replaces  that  obtained  on  the  pre- 
vious sweep.  An  example  of  its  output  is  shown  in  Fig.  2-7  for  the  case  of 


f =5  Hz  Sine  Wave 
m 

Oscilloscope  Data 

Horiz.  = .l  sec/div 
Vert.  = 1 V/div 


Fig . 2 - 7 - Output  of  Signal  Processor  for  FM  Modulated  Input 


an  FM  signal  of  center  frequency  2.0  MHz  (fQ)  modulated  to  +200  kHz  about 
fQ  at  a 5 Hz  rate  sinusoidally. 

A provision  is  included  for  tracking  single  sideband  suppressed  carrier 
signals,  with  an  identification  of  upper  or  lower  sideband  such  that  when  used 
in  conjunction  with  an  acousto -optic  modulator  the  unit  can  discriminate  the 
sign  of  the  Doppler  shift.  The  signal  feedthrough  at  the  translated  frequency 
can  also  be  discriminated  against  digitally,  thus  eliminating  the  need  for  a 
"notch  filter." 

The  raw  spectral  information  (output  of  the  spectrum  analyzer  is  also 
made  available  to  the  Systems  Engineering  Laboratories  (SEL)  810A  data 
logging  minicomputer  which  is  programmed  to  generate  its  own  estimate  of 
the  spectral  peak. 

2.1.4  Data  Recording  and  Display 

Primary:  The  primary  data  gathering  function  is  performed  by  an  SEL 

8 10A  general  purpose  minicomputer.  Data  gathering  by  the  Mobile  Atmos- 
pheric Unit  is  formatted  by  the  computer  software  and  stored  on  magnetic 
tape  for  subsequent  processing  on  the  Univac  1108.  The  SEL  7-track  tape 
control  and  magnetic  tape  units  allow  digital  recording  of  data  at  800  bpi  at 
45  ips,  which  recording  density  is  common  to  the  Univac  1108  I/O  system. 

The  data  logged  by  the  computer  includes; 

• All  scan  volume  location  parameters 

• "Mode  of  operation"  identifier 

• The  instantaneous  line-of-sight  velocity  information 

• The  Doppler  spectrum  peak  strength 

• Full  spectrum  intensity  and  frequency  information  (optical) 

• A data  quality  identifier. 

Properties  of  the  Doppler  spectrum,  namely,  the  amplitude  and  frequency 
corresponding  to  the  spectral  peak  are  obtained  as  a result  of  on-line 


computer  processing  except  for  the  frequency  which  is  also  obtained  by  the 
spectral  peak  locator  (velocity  processor)  discussed  previously.  The  latter 
allows  some  flexibility  for  on-line  operator  displays  (see  below). 


Secondary:  The  velocity  processor  output  estimate  of  the  instantaneous 

line-of-sight  velocity,  updated  at  a 70  Hz  rate,  is  available  in  analog  format 
which  can  be  recorded  directly  on  a strip  chart  recorder,  an  option  which  is 
extremely  useful  during  the  VAD  mode  of  operation  for  monitoring  the  char- 


acteristic profile. 

Some  overall  views  of  the  mobile  unit  hardware  as  utilized  during  this 
program  are  shown  in  Fig s . 2 - 8 through  2-11. 

2.2  WINDS  ALOFT  SENSING 

Using  the  basic  system  outlined  previously  it  is  possible,  by  scanning 
operations,  to  determine  the  three -component  wind  field  at  any  specified 
altitude  between  50  and  2000  feet.  The  scanning  method  employed  is  com- 
monly referred  to  as  the  velocity  azimuth  display  (VAD)  technique  which  was 
first  utilized  by  Lhermitte  and  Atlas  in  conjunction  with  a microwave  radar. 

The  telescope  is  focused  at  the  altitude  of  interest,  the  beam  being  di- 
rected at  a zenith  angle,  (3.  The  beam  is  then  scanned  in  azimuth,  thus  tracing 
out  a circle  at  the  selected  altitude  (Fig.  2- 12). 

The  instantaneous  line-of-sight  component  of  velocity  within  the  sensing 
volume  as  measured  by  the  LDV,  v , is  given  by 

v = v,  sinQ  cos(0  - 0 ) + w cosS 
r h ^ o 

Vj  and  0 respectively,  being  the  speed  and  direction  of  the  horizontal  wind 
motion  and  w the  vertical  motion  at  the  height  being  sampled.  The  azimuthal 

j|c 

Also  known  as  conical  scan  technique  because  of  beam  scanning  configuration. 
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dependence  of  vf  is  sufficient  to  yield  the  horizontal  speed  and  direction  and 
vertical  component  of  velocity,  respectively. 

In  the  present  mode  of  operation,  the  system  is  unable  to  distinguish 
between  positive  and  negative  values  of  v . Therefore,  it  is  the  absolute  value 
°f  vr  ( | v | ) that  is  sensed.  This  results  in  a signal  as  shown  in  Fig.  2-13  in- 
stead of  the  sinusoidal  signal  as  shown  in  Fig.  2- 12.  This  results  in  an  am- 
biguity of  180  deg  in  the  wind  direction  since  it  is  uncertain  which  peak  in 
Fig.  2- 13  represents  looking  into  the  wind.  In  practice,  no  problem  occurs 
because  the  operator  records  approximate  wind  direction,  and  the  data  process- 
ing technique  can  then  calculate  exact  wind  direction.  This  resolves  all  wind 
direction  ambiguities  if  the  operator's  input  estimate  is  within  +89  deg  of  the 
true  wind  direction. 

While  operating  in  the  VAD  mode  the  system  is  capable  of  measuring 
winds  at  n (n  = 1 through  8)  altitudes  (that  can  be  dialed  in  by  using  thumb - 
switches)  in  sequence  over  a total  time  period  of  5np  sec  where 

5 sec  = time  for  conical  sweep  at  one  altitude 

n = number  of  altitudes  to  be  interrogated 

p = number  of  VAD  scans  at  each  altitude 
(can  be  chosen  to  be  1 through  7). 

During  this  investigation  n = 2,  p = 1 were  utilized,  thus  allowing  the  meas- 
urement of  wind  at  two  altitudes  every  10  seconds.  Alternating  between  the 
two  altitudes  (28m  and  43m)  was  performed  for  the  test  duration  required. 


ig.2-12  - Principle  of  VAD  Operati 


Section  3 

DATA  PROCESSING  FOR  WIND  MEASUREMENT 

Acquisition  and  processing  of  the  LDV  signature  is  accomplished  by 
means  of  a compact  data  handling  system  developed  specifically  for  the 
Lockheed-Huntsville  MAU.  The  general  elements  of  the  MALI  data  acquisi- 
tion and  data  processing  system  are  shown  in  Fig.  3-1.  The  digitized  LDV 
intensity  versus  frequency  signal  along  with  its  coordinates  in  space  is  fed 
into  the  SEL  810  minicomputer.  Reprocessing  of  the  LDV  signal  is  carried 
out  on  the  minicomputer  utilizing  on-line  computer  programs  written  in  SEL 
machine  language.  Information  from  the  SEL  810  is  stored  on  magnetic  tape 
and  is  used  as  an  input  to  the  off-line  processing  algorithms.  Off-line  process- 
ing of  the  LDV  signal  is  carried  out  on  a Univac  1108  computer  with  programs 
written  in  FORTRAN  language  and  using  card  inputs  with  information  from  the 
logs  to  supplement  the  data.  The  flow  of  data  from  the  MAU  is  sketched  in 
Fig.  3-2  showing  both  the  on-line  and  off-line  data  processing  routines.  On- 
line manipulation  of  the  data  is  carried  out  by  the  SEL  Data  Logger  program. 
The  off-line  processing  is  carried  out  by  the  VAD  and  Vortex  Track  program. 
The  final  output  consists  of  printouts  and  plots.  A description  of  the  data 
logger,  and  the  VAD  program  and  their  operational  characteristics  is  given 
in  the  following  sections. 

3.1  DESCRIPTION  OF  LDV  SOFTWARE  SYSTEM 

Data  acquisition  in  the  Mobile  Atmospheric  Unit  is  carried  out  by  the 
SEL  Data  Logger  program.  A sweeping  spectrum  analyzer  is  used  to  detect 
the  Doppler  shift  frequency.  A diagram  of  the  output  of  the  spectrum  analyzer 
is  shown  in  Fig.  3-3.  The  output  of  the  spectrum  analyzer  is  the  value  of  signal 
intensity  for  each  of  one  hundred  frequency  bands  spanning  the  entire  frequency 
scale.  The  data  logger  on  the  SEL  computer  records  the  signal  intensify  for 
each  frequency  band  for  which  the  Doppler  frequency  shift  exceeds  the  velocity 
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Fig.  3-3  - Typical  Spectrum  Analyzer  Output  in  VAD  Scan 
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threshold  and  the  signal  intensity  exceeds  the  amplitude  threshold.  This  data 
is  stored  on  magnetic  tape  for  off-line  processing.  A flow  diagram  of  the 
data  logger  is  shown  in  Fig.  3-4. 

An  off-line  conversion  program  identifies  the  frequency  band  with  the 
maximum  amplitude  LD V signal.  This  program  saves  the  maximum  ampli- 
tude l.DV  signal,  Img,  and  its  corresponding  frequency,  vmg.  which  is  above 

both  the  amplitude  and  frequency  thresholds.  The  definition  of  I and  V 

ms  ms 

and  the  shape  of  the  characteristic  LDV  spectrum  is  shown  in  Fig.  3-3.  It 
can  be  seen  that  V is  the  characteristic  velocity  associated  with  the  flow 
phenomenon.  The  output  from  the  off-line  conversion  program  consists  of 
Vms  as  a function  of  time  and  spatial  coordinates  of  the  sensing  volume. 

From  the  output  of  this  program  the  wake  vortex  velocity  field  or  wind  field 
can  be  reconstructed  using  an  additional  off-line  processing  routine. 

The  two-step  process  used  to  generate  ^ as  a function  of  time  and 
space  was  necessitated  by  the  adjunct  purpose  of  relating  laser  attenuation  to 
visibility.  The  output  tape  of  the  data  logger  is  to  be  used  for  the  data  analy- 
sis of  that  study.  In  the  more  usual  situation  of  wind  measurement  (i.e.,  when 
laser  attenuation  data  are  not  needed)  the  two  steps  are  accomplished  in  a one- 
step  process  using  the  SEL  810  computer.  The  output  of  the  SEE  computer 
in  that  case  is  identical  with  that  of  the  off-line  conversion  program  described 
above.  The  flow  chart  for  this  data  logger  is  shown  in  Fig.  3-5. 

Final  processing  of  the  LDV  measurements  is  carried  out  by  the  VAD 

and  Vortex  Track  program.  A macro  flow  chart  of  the  VAD  and  Vortex  Track 

program  is  shown  in  Fig.  3-6.  The  VAD  processing  for  wind  measurement 

and  vortex  tracking  algorithm  are  contained  in  the  same  program.  Vortex 

tracking  was  not  used  in  this  test.  In  this  off-line  program  the  array  of  V 

7 max 

values  which  is  a function  of  time  and  space  is  processed  to  yield  the  three- 
dimensional  wind  field  (VAD  mode)  or  the  aircraft  wake  vortex  trajectories 
(vortex  mode).  The  processing  of  the  VAD  measurements  involves  the  com- 
putation of  the  u,  v,  w wind  components  from  the  characteristic  sinusoidal 


VAD  LDV  signature  discussed  earlier  in  Section  2 and  is  described  in  more 
detail  in  Ref.  5.  The  program  is  geared  to  handle  both  the  translated  and 
non -translated  LDV  signal.  A translated  signal  is  provided  when  the  LDV 
system  includes  a frequency  translator,  which  distinguishes  between  positive 
and  negative  values  of  line-of- sight  velocity.  A non-translated  signal  pro- 
vides only  the  absolute  value  of  line-of- sight  velocity  as  shown  in  Fig. 2-13. 
However,  during  the  course  of  this  research  effort  all  of  the  data  acquisition 
and  data  processing  was  done  in  the  non -tran slate  mode.  Three  basic  tech- 
niques have  been  implemented  to  compute  the  three-dimensional  wind  com- 
ponents as  follows:  (1)  a peak  algorithm  where  the  magnitude  and  location  of 
the  peak  signal  in  the  sinusoidal  LDV  VAD  signature  is  used  to  compute  the 
velocity  components;  (2)  a spectral  processing  for  the  winds  using  the  derecti- 
fied  signal;  and  (3)  a sine  curve  fit.  The  final  output  is  a plot  (and  printout) 
of  the  u,  v and  w velocity  components  as  a function  of  altitude  and  time. 

3.2  DATA  PROCESSING  ALGORITHMS  FOR  WIND  MEASUREMENT 

Lockheed-Huntsville  has  developed  three  algorithms  for  calculating  the 
mean  wind  speed  and  direction  from  the  VAD  signature.  For  each  of  these 
algorithms,  mean  wind  and  direction  is  calculated  for  each  5 second  VAD 
sweep.  Standard  deviations  of  wind  speed  and  direction  can  be  calculated 
from  multiple  VAD  sweeps.  The  data  output  of  the  LDV  system  operating  in 
the  VAD  mode  are  line-of-sight  velocities  measured  at  a selected  number 
(usually  350  in  the  current  Lockheed  system)  of  distinct  points  around  the 
VAD  cone.  We  recall  that  the  line-of-sight  velocity  signature  is  theoretically 
sinusoidal  in  the  VAD  mode  (cf.,  Fig.  2- 12). 

For  all  of  the  processing  algorithms,  preprocessing  of  the  data  occurs; 
this  preprocessing  includes: 

1.  Save  line-of-sight  velocities  for  one  rotation  of  scanner 

2.  If  two  or  more  rotations  occur  at  the  same  altitude,  average  with 
previous  rotations 

3.  Assign  azimuth  angle  to  each  point  (assuming  constant  rotation 
rate) 
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4.  Edit  points  to  eliminate  spurious  points 

5.  Plot  line-of-sight  velocity  versus  azimuth  angle. 

The  edit  criterion  for  the  elimination  of  spurious  points  is  that  the  i*^1 
point  is  eliminated  if 


v ,-v  ...  >.2v 

r,  1 r, i+I  r , i + I 


and 
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A sample  plot  of  unedited  line-of-sight  velocity  versus  azimuth  angle  is  shown 
in  Fig.  3-7. 


3.2.1  Peak  Algorithm 

For  the  calculation  of  wind  velocity  by  the  peak  algorithm,  the  calcula- 
tion procedure  is 


1.  Filter  data  with  an  n point  moving  average 

2.  Identify  the  two  peak  velocity  points 
minimum  of  90  deg  apart 


V . and  V that  occur  at  a 
pi  p2 


3.  Compute  horizontal  component  of  wind  velocity 


V, 


V , + v , 
pi  Pz 
2 sinp 


4.  Compute  horizontal  wind  angle  with  help  of  estimated  wind 
direction 

5.  Compute  vertical  component  of  wind  velocity 


V 


Pi  " P2 
2 cosP 


6.  Derectify  VAD  signal  if  no  translator  is  present  and  plot  derecti  fied 
signal 
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The  individual  data  points  are  filtered  with  an  n point  moving  average 
filter  (n  is  usually  21).  Thus,  each  line-of-sight  point  is  filtered  by 


V 

r,  i 


V .+...+  V . , + V . + V V . /n 
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where  ^ is  the  filtered  value  of  line-of-sight  velocity  to  be  used  in  further 
calculations.  A plot  of  the  filtered  line-of-sight  velocities  for  a 21  point  filter 
is  shown  in  Fig.  3-8.  Additional  samples  of  the  LDV  signature  (including  raw 
data,  filtered  data,  and  derectified  data)  are  presented  in  Appendix  A. 

When  the  LDV  data  are  measured,  an  approximation  of  the  wind  angle  is 
recorded.  The  calculated  wind  angle  is  the  azimuth  angle  of  the  peak  which  is 
closer  to  the  estimated  wind  angle.  The  wind  angle  plus  90  deg  is  the  angle  at 
which  the  line-of-sight  velocity  is  theoretically  zero.  This  angle  is  used  for 
the  derectification  of  the  line-of-sight  signal.  A plot  of  the  derectified  (edited 
but  unfiltered)  line-of-sight  velocity  is  presented  in  Fig.  3-9. 


3.2.2  Fourier  Coefficient  Algorithm 


The  Fourier  coefficient  algorithm  (or  spectral  algorithm)  computes  the 
fundamental  harmonic  of  the  line-of-sight  velocity.  The  Fourier  series  for  a 
generalized  periodic  function  is 
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where  the  Fourier  coefficients  are  given  by 
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whore  the  L discrete  points  are  uniformly  spaced  between  0 and  2tt  and  V. 
is  the  line-of-sight  velocity  for  the  it!l  point  in  the  scan.  For  the  calculation 
of  the  three  wind  components  it  is  necessary  to  calculate  only  aQ,  a^,  and  b^. 
However,  the  second  and  third  harmonics  (i.e.,  n = 2 and  n = 3)  are  also  cal- 
culated as  an  indication  of  the  nonuniformity  of  the  wind.  A correction  factor 
must  be  applied  to  compensate  for  the  absence  of  data  points  wiiich  lie  below 
the  velocity  threshold.  Although  the  line-of-sight  signature  in  the  VAD  mode 
is  theoretically  sinusoidal,  the  existence  of  the  velocity  threshold  causes  the 
signature  to  appear  as  shown  in  Fig.  3-10.  For  the  purposes  of  deriving  the 
correction  factor,  the  origin  is  chosen  so 

a 

v = — ~ + b , sin9 
r 2 1 


In  Fourier  analysis  af)/2  is  the  dc  component. 

If  the  vertical  component  of  velocity  is  small  compared  with  the  magnitude 
of  the  velocity  threshold,  the  angle  for  which  the  signature  is  zero  is  Z.7r,/2 
as  shown  in  Fig.  3-10.  The  parameter,  Z,  is  the  ratio  of  zero  points  to  total 
points  in  the  VAD  scan.  Let  a^and  b^  be  the  values  of  the  Fourier  coeffi- 
cients calculated  from  the  line-of-sight  velocity  signature.  The  integral 
form  of  the  equations  for  Fourier  coefficients  is  used.  For  the  theoretical 
velocity  signature, 

2n 

a'  = 1/jr  f (a  /2  +b  sin9)  <10  = a 
o ' J o'  1 o 
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However  for  the  actual  velocity  signature 
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Therefore  a correction  factor  of  l/(l-/,)  must  be  applied  to  the  calculate 
value  of  a . 
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Fig.  3-10  - Line-of -Sight  Velocity  Signature  of  VAD 
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For  the  sine  term,  for  the  theoretical  velocity  signature 


fc>i  = 1/ 7r  J"  (a^/2  +bj  sin9)  sinQ  d9  = 


For  the  actual  velocity  signature 

■n-Zir/  2 2?r-Zff/2 

b j = 1/ 7T  J (ao/2 +b1  sin9)  sin9  d9  + J'  (ao/2  +b  t sin9)  sin9  d9 

.Z*/2  TT+Zff/2 


Therefore,  a correction  factor  l/[l-Z+  (sin27r)/^]  must  be  applied  to  the 
horizontal  wind  components. 

The  calculation  procedure  for  the  spectral  algorithm  is: 

1.  Compute  the  Fourier  coefficients  for  the  first  three  harmonics 
(n  = l,  2,  3) 

L 

a^  = 2/L  V cos(2?rin/L) 

i = 1 

L 

b^  = 2/L  V.  sin(27rin/L) 

i = l 

2.  Compute  corrected  fundamental  harmonic 

Corrected  Value  = - — ^ ^ — r-y—  * Calculated  Fundamental 

I - z.  + (sinz7T)/7r 

3.  Compute  vertical  wind  correction 


Corrected  Value  = -j — p * Calculated  Value 


4.  Compute  horizontal  velocity 

v -S3. 

h sin(3 

5.  Compute  horizontal  angle 

Angle^  - Atan(bj/aj) 

6.  Compute  vertical  wind  velocity 


w 


-a 

o 


2 cosj3 


3.2.3  Sine  Curve  Fit  Technique 

The  line-of- sight  velocity  signature  is  sinusoidal  for  a uniform  wind. 
Therefore,  we  find  the  best  sinusoidal  wave  which  fits  the  data  in  a least 
squares  sense.  Therefore  we  wish  to  minimize 


'y  ^ (V.  - C - A cos9.  - R sinfM^ 


where  V\  is  the  line-of-sight  velocity  (derectified)  at  point  i,  0.  is  azimuth 
at  point  i.  Thus  we  obtain  A,B,  and  C by 
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The  steps  for  calculating  wind  using  the  least  squares  algorithm  are: 

1.  Find  least  squares  curve  fit  for  a sine  wave  to  the  data 

Minimum  = £ (V.  - C - A cosGj  - B sin^)2 

wher  e 

Y.  is  line-of- sight  velocity  at  point 
0.  is  azimuth  at  point 

C,  A and  B are  coefficients  to  be  solved  for. 

2.  Compute  horizontal  velocity 

h sin(cone  angle/2) 

3.  Compute  horizontal  angle 

Angle^  = Atan(B/A) 

4.  Compute  vertical  wind  velocity 


a."— _ 


Section  4 

DATA  COLLECTION  AND  ANALYSIS 

In  this  section  the  data  collection  procedure  and  interpretation  of  the 
data  are  described.  Proper  interpretation  of  the  data  is  very  important. 

The  Lockheed  Mobile  Atmospheric  Unit  was  built  as  an  experimental  unit 
for  the  measurement  of  atmospheric  wind  phenomena.  At  the  time  of  the 
beginning  of  the  fog  tests  at  Otis  AFB,  the  Lockheed  unit  had  been  in  opera- 
tion for  approximately  ten  months.  During  this  time  period,  much  informa- 
tion about  the  operation  of  remote  sensing  systems  was  gained,  and  new 
operating  techniques  were  generated.  For  example,  the  least  squares  sine 
algorithm  was  originated  because  of  inaccuracies  in  the  peak  and  spectral 
algorithms  observed  in  field  tests. 

During  the  field  tests,  design  improvements  in  both  system  hardware 
and  software  have  become  apparent,  and  many  of  these  improvements  have 
been  incorporated  into  LDV  systems  built  for  the  U.S.  Army  and  the  Depart- 
ment of  Transportation.  Some  of  these  improvements  are  discussed  in  this 
section  as  an  indication  of  the  manner  in  which  data  resulting  from  a system 
designed  for  use  in  a fog  dispersal  system  would  differ  from  the  data  pre 
sented  herein . 

The  data  described  herein  is  the  first  data  collected  in  fog  with  the 
laser  Doppler  velocimeter.  As  expected,  the  system  performed  well.  How  - 
ever, certain  operational  and  design  improvements  could  further  enhance  the 
quality  of  the  data  in  a system  designed  for  use  in  fog.  The  improvements 
are  discussed  in  this  section. 

4.1  TEST  DESCRIPTION 

The  Lockheed  Mobile  Atmospheric  Unit  was  operational  at  the  Otis 
AFB  test  site  from  7 September  through  30  September  1 976.  The  unit  was 
placed  on  standby  status,  and  the  unit  was  ope-ated  for  wind  measurement 
during  periods  of  fog.  A diagram  of  the  Otis  AFB  test  site  is  shown  in 
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Fig. 4-1.  For  the  test,  the  direction  of  the  towers  was  assumed  to  be  north. 
The  positive  u coordinate  was  taken  as  wind  from  the  west,  and  the  positive 
v component  was  taken  as  wind  from  the  south.  A table  of  the  days  and  times 
on  which  data  were  taken  is  shown  in  Table  4-1.  During  the  test  runs,  the 
altitude  at  which  the  VAD  scans  were  made  alternated  between  30  m and  45  m. 
Thus,  in  each  minute  time  period,  six  scans  were  made  at  30  m and  six  scans 
were  made  at  45  m.  Each  scan  requires  five  seconds  of  time. 

The  brief  time  breaks  in  the  data  are  for  two  reasons.  The  first  was 
to  adjust  the  laser  power  for  the  adjunct  purpose  of  obtaining  data  on  laser 
attenuation  in  fog  and  was  unrelated  to  the  primary  purpose  of  the  test.  The 
laser  power  fluctuates  slightly  with  time.  In  the  condition  for  which  wind 
data  only  are  required,  laser  power  fluctuation  is  never  large  enough  to  be 
of  concern.  However,  for  laser  attenuation  studies,  it  was  believed  that  the 
power  should  be  kept  constant. 

The  second  reason  for  the  short  interruptions  was  the  removal  of 
moisture  condensation  from  the  mirrors  in  the  scanner.  In  a system  de- 
signed for  operation  in  fog,  a small  heater  in  the  scanner  would  prevent 
such  condensation. 

4.2  DATA  PRESENTATION 

Sample  data  output  are  shown  in  Figs.  4-2  and  4-3.  A complete  tabular 
listing  of  the  data  measured  during  the  test  is  presented  in  Appendix  13.  Hie 
date  is  the  test  site  date  on  which  the  run  was  initiated.  All  times  are 
Greenwich  Mean  Time.  The  height  is  the  indicated  altitude  (meters)  on  the 
laser  focusing  system.  Before  each  test  run,  the  laser  ranging  system  is 
calibrated  to  assure  that  the  actual  altitude  of  the  data  is  the  desired  altitude. 
Thus,  an  indicated  altitude  of  28  m corresponds  to  an  actual  altitude  of  30  m, 
and  an  indicated  altitude  of  43  m corresponds  to  an  actual  altitude  of  45  m. 
The  data  wore  reduced  to  give  five  minute  averages  beginning  at  even  five 
minute  time  periods.  Although  the  run  of  16  September  started  at  09:15;  30, 
the  first  even  five  minute  time  period  began  at  09:20:00.  I lie  following  in- 
terpretations are  placed  on  each  of  the  columns. 
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Tower  3 
(5  m) 


38  True 
(Called  0 
for  test) 


Fig  . 4 - 1 - Diagram  of  Otis  A I' 13  Test  Site 
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Table  4-1 

DATA  COLLECTION  AT  OTIS  AFB 


Date 

Start 

Stop 

Thr  e shold 
Velocity 
(m/ sec ) 

Approx. 
Wind  Speed 
(m/sec ) 

Approx. 

Wind 

Direction 

! 

Fog  Type 

9-16-76 

9: 15:30 

10:00:00 

0.53 

2 

140° 

Radiation 

10:03:15 

10:30:00 

0.53 

1.5 

c 

o 

o 

rj 

10:31:00 

1 1:00:00 

0.53 

1.5 

200° 

1 1:2  7:20 

1 1:36:4  5 

0.53 

1.5 

200° 

11:41:20 

12:00:05 

0.5  3 

2 

ro 

4- 

O 

0 

I 

12:07:00 

12:12:00 

0.53 

9-18-76 

3: 16:20 

6:08:40 

0.79 

3-5 

150° 

Advection 

6:16:00 

6:51:00 

0.79 

5 

O' 

o 

0 

9-18-76 

7:16:00 

8: 12:00 

1 .06 

7 

180° 

1 9-28-76 

00:25:00 

00:39:00 

1.06 

00:44:00 

03:00:00 

1.06 

j 9-28-76 

3:15:17 

3:25:47 

2.12 

3:  30:00 

3:34:45 

2.12 

3:39:00 

4:41:00 

2.12 

1 
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Fig.  4-3  - Sample  Data  Output  at  156  ft  Altitude 


4.2.1  One  Minute  Means 


One  minute  means  were  calculated  for  each  of  the  minutes  of  the  five 
minute  time  period.  The  minutes  of  the  time  period  are  tabulated  in  the 
first  column.  Each  of  the  numbers  in  this  group  is  an  average  of  the  six 
VAD  scans  that  occur  at  the  given  altitude  during  one  minute. 

Peaks 

The  one  minute  mean  u,  v and  w components  of  wind  as  calculated 
by  the  peak  algorithm  are  tabulated.  The  w component  is  the  vertical  com- 
ponent and  is  positive  vertically  upward.  Speeds  are  given  in  cm/sec  and 
angles  are  given  in  degrees. 

The  data  calculated  from  the  basic  VAD  data  are  horizontal  wind  speed, 
wind  direction,  and  vertical  component  of  wind  velocity.  The  u and  v com- 
ponents are  calculated  from  horizontal  wind  speed  and  direction.  The  speed 
is  the  one  minute  average  of  the  horizontal  component  of  velocity,  and  "TH" 

is  0 the  direction  of  the  horizontal  component  of  wind.  The  u and  v aver- 
o’ 

ages  are  calculated  by  averaging  the  u and  v values  of  the  individual 
VAD  scans  rather  than  from  averaged  values  of  u and  v.  Thus 

u.  = -S.  sinQ 

1 1 0,1 

and 

v.  = -S.  cos9 

1 1 o,  1 

where  the  i subscript  refers  to  the  individual  VAD  scans.  The  average  for 
n scans  in  the  one  minute  average  is 

S = £ S./n 


| ,-c  use  of  :i  procc-sin.-  err  r (since  corrected)  the  w component  should 
be  multiplied  by  -O.B. 


L 


and 


It  is  noted  in  general  that 


v = 


S0o,/n 

£u/n 

Evn 


u / -S  sin9 
' o 

V / - S COS0 


and 


e , J-Z  , -2 

S p * u + v 


The  same  averaging  characteristic  results  from  resolving  cup  and  vane 
anemometer  data  into  components  and  averaging  or  when  comparing  averaged 
data  from  cup  and  vane  anemometers  with  averaged  data  from  propeller  ane- 
mometers. The  characteristic  is  particularly  pronounced  in  light  and  variabli 
wind  conditions. 

Fourier  Coefficients 

Results  are  shown  for  wind  calculations  performed  using  the  spectral 
algorithm.  The  symbols  have  the  same  meaning  as  for  the  peaks  algorithm. 
The  "2D"  is  the  magnitude  of  the  second  harmonic  and  is  the  average  of 

VZ  2 

a , + b^  where  a^  is  the  second  Fourier  cosine  coefficient  and  \>^  is  the 
second  Fourier  sine  coefficient.  Similarly,  "3D"  is  the  average  of  + b*'  . 

The  magnitude  of  these  numbers  is  an  indication  of  turbulence  with  the  fre- 
quency of  the  turbulence  increasing  with  the  number  of  tin-  harmonic. 

Sine  Wave  Fit 


Results  are  shown  for  the  least  squares  sine  fit  algorithm  with  similar 
meanings  of  the  symbols  as  for  the  other  two  algorithms.  The  "SP  is  an 
indication  of  the  deviation  of  the  data  points  from  a perfect  sine  fit. 


Because  of  a processing  error  (since  corrected)  the  w component  should 
lie  multiplied  by  -1. 

■ 


w 


4.2.2  Cumulative  Means 

The  cumulative  means  are  averages  of  individual  VAD  scans  taken 
over  the  time  period  indicated.  For  example,  the  three  minute  cumulative 
mean  is  the  average  of  the  VAD  scans  taken  over  the  first  three  minutes  of 
the  five  minute  time  period.  They  are  not  averages  of  the  one  minute  means. 
The  two  averaging  methods  would  be  identical  if  the  number  of  VAD  scans  in 
each  one  minute  average  were  always  identical.  The  five  minute  cumulative 
means  are  the  averages  over  the  five  minute  time  period  of  interest. 

4.2.3  One  Minute  Standard  Deviations 

This  section  lists  the  standard  deviation  of  each  of  the  variables  for 
each  one  minute  time  period.  Each  standard  deviation  is  the  standard  de- 
viation of  the  six  VAD  scans  in  the  one  minute  period. 

4.2.4  Cumulative  Standard  Deviations 

The  cumulative  standard  deviations  are  analogous  to  the  standard 
means.  They  are  the  standard  deviations  of  the  individual  VAD  scans  taken 
over  the  time  period  indicated.  In  general,  the  cumulative  averages  of  the 
standard  deviations  are  not  averages  of  the  one  minute  means  because  of 
the  mathematical  definition  of  standard  deviation. 

4.3  DATA  ANALYSIS 

Although  meteorological  tower  data  are  not  available  for  direct  com- 
parison with  laser  measured  winds,  a discussion  of  previous  wind  comparisons 
with  meteorological  towers  is  helpful  in  the  interpretation  of  the  results  of 
this  test.  Observations  about  operation  in  fog  are  also  appropriate. 

4.3.1  Comparison  of  Computational  Algorithms 

The  relative  accuracies  of  the  three  computational  algorithms  can  be 
seen  by  comparison  with  meteorological  flat  a in  past  tests. 
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The  winds  calculated  with  the  three  algorithms  described  in  Section  3 
are  shown  in  Figs.  4-4  through  4-6  for  data  measured  at  the  Wave  Propaga- 
tion Laboratory,  National  Oceanic  and  Atmospheric  Administration,  Boulder, 
Colorado.  The  winds  were  measured  in  March  1976.  The  sole  difference 
between  the  three  parts  of  each  figure  is  the  data  processing  algorithm  used. 
Some  general  observations  are  that  the  peak  algorithm  tends  to  give  mean 
wind  values  that  are  biased  slightly  high  compared  to  the  tower  measured 
mean  winds.  This  is  primarily  due  to  temporal  variations  in  the  wind  at 
higher  frequencies  than  the  fundamental  frequency  in  the  VAD  signature 
thus  the  peak  of  the  signature  (i.e.,  the  fundamental  frequency  with  super- 
imposed higher  frequency  variations)  has  a value  that  is  higher  than  the 
averaged  values.  This  interpretation  is  confirmed  by  the  runs  for  which  the 
greatest  deviation  above  the  45  degree  line  occur  are  the  runs  for  which  the 
greatest  value  of  wind  standard  deviation  (measured  by  the  tower)  occur.  In 
particular,  for  the  points  of  Run  9 farthest  off  the  45  degree  line  in  Figs. 4-5 
and  4-6,  the  standard  deviation  of  wind  speed  is  1.82  m/sec  and  1.75  m sec, 
respectively,  whereas  the  standard  deviation  of  wind  speed  was  usually  be- 
tween 0.9  m/sec  and  1.4  m/sec. 

It  is  also  observed  that  the  spectral  algorithm  tends  to  give  mean  wind 
values  that  are  biased  slightly  low  compared  to  the  tower  measured  mean 
winds.  It  was  observed  during  the  data  processing  of  the  Boulder  data  that 
the  magnitude  of  the  bias  is  a function  of  the  number  of  valid  line  of  si.dit 
velocities  obtained  during  the  VAD  scan.  The  bias  is  in  significant  for  data 
where  more  than  50  valid  data  points  are  obtained  during  each  YA1)  scan,  t 
is  noted  that  a number  of  data  points  below  this  value  was  obtained  only  rarely 
at  Table  Mountain  and  that  the  low  number  of  data  points  would  never  occur  m 
air  with  a larger  quantity  of  natural  aerosol  than  the  very  clean  air  of  Table 
Mountain.  Tire  least  squares  sine  algorithm  seems  to  present  no  systematic 
bias,  and  the  scatter  of  the  data  about  the  45  degree  line  is  less  >r  tin-  sine 
algorithm  than  for  the  other  two  processing  techniques.  A comparison  of  tin 
peak  algorithm  and  the  spectral  algorithm  with  the  sine  algorithm  for  the  data 
measured  at  Otis  AFB  is  shown  in  Fig.  4-7.  In  similar  figures  in  the 
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W ind  Speed  from  Sine  Algorithm  (m/sec) 


7 - Comparison  of  Datii  Processing  Algorithms  Measured 
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Boulder  data  (Figs.  4-8  and  4-9)  it  was  observed  that  the  speed  calcu 
lated  by  the  peak  algorithm  was  almost  always  greater  than  that  calculated 
by  the  sine  algorithm.  The  same  observation  is  valid  for  the  data  greater 


than  4 m/sec  at  Otis  AFB,  but  for  light  winds  (<  2 m/sec)  the  peak  algorithm 
often  gives  a value  which  is  less  than  that  derived  from  the  sine  algorithm. 

4.3.Z  Variations  in  Fog  Density 

One  of  the  problems  associated  with  wind  measurement  in  fog  is  meas- 
urement in  both  fog  and  clear  air.  Unfortunately,  the  magnitude  of  this  prob- 
lem was  not  completely  anticipated  and  some  data  were  lost.  In  fog,  for  the 
manner  in  which  the  data  were  taken,  the  backscatter ed  signal  is  very  large 
and  must  be  attenuated  to  avoid  going  off  scale  in  the  spectrum  analyzer. 

With  this  quantity  of  attenuation,  when  the  focal  volume  hits  clear  air.  an  ade- 
quate signal  cannot  be  obtained.  This  problem  occurs  if  the  data  are  taken 
on  a linear  scale.  The  Lockheed  system  also  has  a logarithmic  scale  which 
avoids  the  problem  by  displaying  and  recording  the  logarithm  of  the  intensity 
of  the  return  signal.  Since  the  magnitude  of  the  problem  was  not  anticipated, 
all  data  during  the  test  were  taken  in  the  linear  mode. 

The  problem  is  illustrated  by  the  data  for  18  September.  From  the 
magnitude  of  the  wind  it  is  apparent  that  the  fog  is  an  advection  fog.  Valid 
data  for  wind  speed  on  the  order  of  3 m/sec  to  4 m/sec  is  obtained  at  both 
altitudes  from  3:20  to  5:20.  However,  after  5:20  the  wind  increased  to  an 
order  of  5 m/sec  to  6 m/sec.  This  would  cause  the  advection  fog  to  lift 
slightly,  providing  clearer  air  at  the  lower  altitude.  Tin-  data  taken  at  30  m 
altitude  between  5:20  and  6:50  is  therefore  somewhat  sparse.  During  the 
brief  interruption  from  6:51  to  7:16  the  laser  power  and  signal  attenuation 
were  adjusted,  and  valid  data  were  obtained  thereafter. 

A similar  condition  is  seen  in  the  data  of  16  September.  The  wind 
speed  indicates  a radiation  fog.  Valid  data  are  obtained  at  the  beginning  of 
the  run.  As  the  depth  of  the  fog  decreased  after  sunrise,  the  system  began 
to  lose  data  at  the  higher  altitude.  Adjustment  of  laser  power  and  attenuation 
at  11:40  restored  a good  signal  at  both  ^lltitudes. 
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Fig.  4 -8  - Comparison  of  Data  Processing  Algorithms  at  60  m Altitude 
Measured  at  Table  Mountain,  Colorado 


Hg.4-Q  - Comparison  of  Data  Processing  Algorithms  at  150  m Altitud 
Measured  at  Table  Mountain,  Colorado 
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As  mentioned  previously,  these  adjustments  would  have  oeen  unneces- 
sary if  the  laser  intensity  had  been  measured  in  the  logarithmic  scale.  In 
order  to  verify  that  an  adequate  signal  could  be  obtained  both  in  and  out  of 
fog  conditions,  Lockheed-Huntsville  made  measurements  in  low  cloud  condi- 
tions at  its  Huntsville  facility  on  7 December  1976.  The  reported  meteoro- 
logical conditions  at  the  Huntsville-Madi  son  County  Jetplex  (10  miles  from 
the  test  site)  are  shown  in  Table  4-2,  and  wind  conditions  are  shown  in  Figs. 
4-10  through  4-17.  More  complete  data  from  this  set  of  measurements  are 
shown  in  Appendix  C.  Adequate  data  were  obtained  for  both  in  and  out  of  the 
clouds  using  the  logarithmic  scale. 

4.3.3  Light  and  Variable  Winds 

The  existence  of  light  and  variable  winds  presents  some  unique  condi- 
tions which  must  be  considered  in  the  data  analysis.  The  first  of  these  is 
important  when  comparing  the  LDV  data  with  propeller  anemometer  data. 

For  the  propeller  anemometer  data,  the  wind  can  fall  below  the  threshold 
value  in  one  of  the  component  directions.  However,  because  the  LDV  oper- 
ates in  speed-direction  coordinates  and  then  transforms  the  data  into  ortho- 
gonal coordinates,  it  has  no  threshold  value  in  the  orthogonal  directions. 

Thus  the  LDV  can  indicate  values  of  wind  which  are  below  the  threshold  value 
of  the  anemometers.  This  can  affect  the  averaging  over  a period  of  time. 

In  the  mode  of  LDV  operation  without  a frequency  translator,  the  system 
is  unable  to  distinguish  between  positive  and  negative  velocities.  Therefore, 
at  the  beginning  of  each  run,  the  operator  makes  an  estimate  of  wind  direction, 
and  the  direction  is  chosen  as  the  angle  of  the  peak  of  the  VAD  scan  (cl.  tigs. 
3-7  and  3-8)  which  lies  closest  to  the  estimated  wind  direction.  In  practice 
this  causes  no  problem  if  the  actual  wind  direction  is  within  H49  deg  of  the 
estimated  wind.  In  all  previous  applications  of  the  VAD  technique,  this 
criterion  has  been  easily  met,  and  no  ambiguity  has  resulted  from  the  sys- 
tem limitation.  However,  in  light  and  variable  winds,  the  variation  in  wind 
direction  from  the  estimated  wind  direction  can  occasionally  exceed  90  deg 
for  one  or  more  scans  in  the  one  minute  average.  For  each  scan  for  which 
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Fig.  4- 10  - Wind  Data  Measured  at  30m  Altitude  with  Low  Clouds 
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ind  Data  Measured  at  270  m Altitude  with  Low  Clouds 


Table  4-2 


METEOROLOGICAL  CONDITIONS  FOR  WIND  MEASUREMENT 
IN  LOW  CLOUDS,  HUNTSVILLE -MAD  ISON  COUNTY  AIRPORT 


8:55  CST 

Clouds:  800  ft  Thin  Broken,  1600  ft  Overcast 
Visibility:  2-1/2  Miles  in  Light  Drizzle  and  Fog 
Wind:  320  deg  at  10  knots 


9:55  CST 

Clouds:  400  ft  Thin  Broken,  800  ft  Overcast 
Tops  Reported  at  7000  MSL 
Visibility:  1-1/2  Miles  in  Light  Drizzle  and  Fog 

Wind:  330  deg  at  14  knots 


10:55  CST 

Clouds:  600  ft  Overcast 

Tops  Reported  at  7200  MSL 
Visibility:  2 Miles  in  Light  Drizzle  and  Fog 
Wind;  340  deg  at  12  knots 


this  condition  occurs,  the  wind  angle  is  erroneous  by  exactly  180  deg.  The 
wind  speed  is  correct.  The  u and  v components  are  correct  in  magnitude 
for  each  VAD  scan;  but  are  wrong  in  sign. 

The  one  minute  averages  for  which  this  condition  occurs  can  be  clearly 
identified  in  the  data.  Since  it  is  known  that  there  are  six  scans  in  every 
one  minute  average  and  the  angle  of  one  or  more  of  the  scans  is  wrong  by 
180  deg,  the  standard  deviation  of  the  one  minute  mean  of  the  wind  angle  will 
exceed  60  deg  if  the  condition  occurs.  Similarly,  the  standard  deviation  of 
both  the  u and  v components  must  exceed  74%  of  their  respective  mean 
values.  To  see  this,  assume  that  the  true  value  from  each  scan  of  a one 
minute  average  is  a.  Therefore,  for  6 scans  in  the  one  minute  average, 
the  expected  value  and  expected  value  of  the  square  are 

E(X2)  = 6a2/6  E(X)  = 6a/6 

and  the  Variance  is 

2 2 2 2 
E(X  ) - E (X)  = a - a = 0 

If  one  of  the  six  values  has  the  wrong  sign, 

E(X  2 ) = 6 a 2/6  E(X  ) = 4a/ 6 

and  the  variance  is 

E(X2)  - E2(X)  = a2  - (4a/6)2  = 0.52a2 

or  the  standard  deviation  is  0.74  a.  Thus  the  condition  described  above  is 
suspected  if  the  standard  deviati  ons  of  both  the  u and  v components  exceeds 
74%  of  their  respective  mean  values  f . » r tin-  one  minute  average. 


d 


In  some  conditions  (e.g.,  43  m data  of  14  September  1976  ending  at 
10:40:00)  the  standard  deviation  for  winds  c£ilculated  from  the  sine  algorithm 
may  indicate  that  this  condition  exists,  whereas  the  standard  deviation  for  the 
peak  algorithm  indicates  that  the  problem  does  not  exist.  When  this  condi- 
tion occurs,  the  wind  data  calculated  with  the  peak  algorithm  can  be  expected 
to  give  more  reliable  data  than  the  wind  data  calculated  with  the  sine  algo- 
rithm. 

It  is  noted  that  this  condition  occurs  only  rarely  for  the  light  winds  of 
16  September  and  never  for  the  higher  winds  of  18  September.  The  higher 
winds  are  less  susceptible  to  wind  speed  changes  of  more  than  90  deg. 

This  limitation  of  the  Lockheed  LDV  as  configured  for  the  test  has  been 
recognized  in  the  past.  A translator  is  currently  being  installed  in  the 
Lockheed  system,  and  translators  have  been  installed  in  the  LDV  systems 
built  for  the  U.S.  Army  Missile  Command  and  Department  of  Transportation. 
The  presence  of  a translator  removes  all  ambiguity  associated  with  the  di- 
rection of  the  wind  and  solves  the  problem  described.  It  also  obviates  the 
need  for  estimates  of  wind  direction  to  be  made  by  the  system  operator.  The 
absence  of  a translator  has  caused  no  problem  in  past  measurements  of  wind 
using  the  VAD  mode.  However,  it  is  clear  that  a translator  will  be  necessary 
for  measurement  of  wind  in  light  and  variable  conditions. 
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Section  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  ability  of  the  laser  Doppler  velocimeter  to  measure  winds  in  fog 
has  been  demonstrated.  The  continuity  of  the  data  recorded  and  the  quantity 
of  data  recorded  were  adversely  affected  by  the  secondary  purpose  of  the 
test.  This  purpose  was  the  measurement  of  laser  attenuation  in  fog  and  re- 
quired that  the  laser  power  be  kept  constant  within  limits  more  restrictive 
than  required  by  wind  measurements  alone. 


The  addition  of  a translator  (to  distinguish  between  positive  and  negative 
line -of -sight  velocities)  and  real  time  data  processing  will  further  enhance  the 
ability  of  the  LDV  to  measure  winds  in  fog.  Both  of  these  additions  are  planned 
lor  1977.  The  ability  of  the  LDV  to  measure  wind  in  patchy  fog  or  in  fog  and 
clear  air  above  the  fog  is  enhanced  by  the  use  of  the  logarithmic  scale  of  laser 
intensity  rather  than  the  linear  scale. 


It  is  recommended  that  the  future  activities  of  the  Lockheed  LDV  and 
the  Department  of  Transportation  LDV  be  monitored.  It  is  likely  that  both 
of  these  units  will  take  data  in  fog  in  their  normal  use. 
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Appendix  B 

TABULAR  DATA  FOR  WIND  MEASURED  AT  OTIS  A F B , 
MASSACHUSETTS,  IN  FOG  CONDITIONS 
DURING  SEPTEMBER  1976 
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TABULAR  DATA  FOR  WIND  MEASURED  IN  LOW 
CLOUD  CONDITIONS  AT  LOCKHEED  MISSILES 
& SPACE  COMPANY,  HUNTSVILLE, 
ALABAMA,  ON  7 DECEMBER  1976 
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